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'^ . Abstract: The best 3-channel projectile-inelastic [H(ls)+Ps(ls,2s,2p)] 

^ I close-coupling approximation (CCA) is performed for positronium 

I (Ps) and hydrogen (H) collision considering both the atoms in 

r-| ' ground states at the incident channel. The s-wave elastic phase 

0-t- shifts and cross sections in the singlet and triplet channels and 

i ! the total (or integrated) cross section are studied below inelas- 

Q I tic threshold. Resonances in singlet channel using 2-channel and 

c^ I 3-channel projectile-inelastic CCAs confirm earlier prediction [1]. 

c/3 ■ A resonance in triplet channel using 3-channel projectile-inelastic 

.y . CCA is a new addition. 

^ : PACS No.s: 36.10.Dr, 01.55.+b, 03.65.Nk, 34. 

^ ; 1 Introduction 

\0 . To learn new physics, the Ps-H is an ideal system since the target (H) and 

the projectile(Ps) both are hydrogen-like systems with a single electron; 

their wavefunctions are known exactly. Being a bound system of a par- 

T-H ' ticle(electron) and its antiparticle (positron), the charge and mass centers 

^ ■ coincide in positronium(Ps); this gives rise to a zero mean static potential. 

;• . In addition the eight times higher polarizability than hydrogen(H) makes Ps 

.^ ! an interesting probe for new investigations. The theoretical studies of such 

r> I a system is rigorous due to the presence of four charge centers. At very low 

cd ' energies near to zero although elastic scattering is the only active channel 

but the excitation and ionization thresholds are very close. So the quantum 

mechanical influence of different inelastic channels are important and should 

be taken into consideration. It is a challenging problem to get accurate low 

energy data for such a system. 

In collision physics, resonance is an important phenomenon. When a 
microscopic moving object that is a wave, enters into the scattering chamber 
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near the target, it faces interactions. When it comes out of the scattering 
zone, the original incident wave gathers a phase shift and the new wave is 
known as scattered wave. The change in phase which is named as phase 
shift, is the parameter that carries the information of the scattering process. 
A rapid change in phase shift by tt radian in a very narrow energy interval of 
the incident wave is an indication of the existence of a resonance. It carries 
the information of a bound system if in the s-wave elastic scattering and 
below threshold of excitation. One can calculate the width of resonance to 
get the life time (r) of the newly formed system. 

It is an extremely difficult job to detect a resonance since successful iden- 
tification needs (i) a very accurate calculation and (ii) sufficient computation 
facilities. A large number of mess-points in a very small energy interval, gen- 
erally ~ 10^^ — 10"'^ eV is required. It necessitates a high-speed computer 
with a sufficient memory, the knowledge of mathematical computation and 
programming languages e.g. FORTRAN. 

We are interested to discuss the scattering processes at low energies. At 
high energies the projectile usually ignores all the important delicate interac- 
tions with target, so generally it carries no valuable information. At very low 
energies below excitation threshold, elastic scattering is the only real process. 
But from just above the threshold different excitation, ionization etc. chan- 
nels start to open. It is a very difficult task to study the scattering processes 
at intermediate energies due to presence of many different channels and a 
very close coupling among them. The total cross section is the sum of the 
integrated cross sections of all these channels. Besides these, the partial wave 
contribution from higher angular momenta start to dominate since above the 
threshold. Many close peaks may arise in the cross section at the interme- 
diate energies due to opening of different channels and their mixing. These 
are also named as resonances and studied by Higgins & Burke [2-3], Sarkar 
et al [4-5]. However Zhou et al [6] commented against such resonances. So 
the subject of above threshold resonances needs more investigation. Above 
threshold, the phase shift becomes a complex quantity and known as eigen 
phase shift. 

1 . 1 Resonances 

The kind of resonances we are interested in, is feasible only below inelastic 
threshold. At this energy region, the integrated elastic cross section {a^i = 
af^Q + o-f^i + (t£L2 + •••) is the total cross section {(Trotai = (^ei + J2 (^inelastic) 



due to presence of only elastic channel; again s-wave (i.e. L = partial wave) 
dominates. So effectively, the projectile angular momentum L = partial 
wave carries almost all the information of the scattering process. 
The phase shift 5i can be decomposed as 

^i = ii + Vi- 
^i corresponds to the hard sphere scattering or non-resonant part; it does not 
depend on the shape and depth of the potential. The term rji depends on the 
details of the potential. The quantities ^i and r]i vary in general slowly and 
smoothly with the incident particle energy. But in certain cases r]i may vary 
rapidly in a small energy interval of width T about a given energy value Er 
such that we can write 

r 
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In that energy interval the phase shift is therefore given approximately by 

The physical significance of a narrow resonance can be inferred by exam- 
ining the amplitude of the radial wave function inside the interaction region. 
The probability of finding the scattered particle within the potential is much 
larger near the resonance energy E = Eji, so that in that case the particle 
is nearly bound in the well. Thus the resonance may be considered as a 
metastable state whose lifetime r, which is much longer than a typical colli- 
sion time, can be related to the resonance width F by using the uncertainty 
relation AtAE > h. Thus, with At ~ r and AE ~ F, we have r ~ ^. 

The shape of the cross section curve near a resonance as a function of 
energy depends on the non-resonant phase shift C,i [7,8]. For the s-wave 
scattering it is 
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sin%{ER - Ef + cos^iiT^ + sin2ii{ER - E)T /2 
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Two limiting cases for non-resonant phase shift are and 7r/2. In the 
first case the above equation becomes 
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which is symmetric and represents a rise in cross section at the resonance 
energy. In the other case 



^' (i?« -EY + ^i 

which is also symmetric but goes down to zero at the resonance energy. If 
the non-resonant phase shift gets some other value then all sorts of forms of 
the cross section can occur. 

However one needs to investigate both the phase shift and the cross sec- 
tion, but should be careful to use the formulation which derives the cross 
section directly from amplitude and not from phase shift. If resonance is 
true it should be reflected both in phase shift and in cross section. It is 
actually to make a precautionery measure for successful detection since the 
calculation of phase shift involves a tan-inverse function which may create a 
numerical error. 

1.2 Close coupling approximation (CCA) 

The CCA is a highly successful theory to study low energy scattering phe- 
nomenon in atomic physics. The theory was predicted by Massey [9] and 
applied in e^ - atom scattering. Later Burke et al [10] successfully used this 
theory for e^ - atom scattering. Now a days many different groups of the 
world are using this theory for e~ je^ - atom scattering. It was Fraser [11] 
who used it first for Ps - H scattering. Ray & Ghosh [12-13] merits the 
credit because they are the first who supplied detailed and converged results. 
They used a momentum space formalism introduced by Calcutta group [14] 
whereas Fraser used a coordinate space formalism to write the coupled in- 
tegral equations. They again add more channels in the CCA basis [15-28]. 
The studies of Fraser [11,29-30], Fraser et al [31], Hara et al [32] were con- 
fined to static-exchange model i.e. considering only the elastic channel in the 
basis and the H and He targets. Ray [16-18] extended the CCA theory in 
Ps and lithium (Li) scattering using the static-exchange and a two-channel 
CCA models. 

The theory is based on the very basic principle of quantum mechanics 
i.e. the eigen state expansion (ESE) methodology in which the total wave 
function of a quantum mechanical system is expressed as a linear combination 
of all possible states known as basis set. So one has to use a wide channel 
space, but practically it is impossible. The number of unknowns exceed 
the number of equations when non-spherical orbitals like p- and d- states 



of an atom are considered in channel space. So arised the necessity of an 
approximation. We should conserve the total angular momentum quantum 
numbers e.g. 'J' and 'M'. It makes the equations closed i.e. the number of 
unknowns are equal to the number of equations. The accuracy of the method 
depends on the choice of basis set. It should be remembered that when 
we study a particular channel (e.g. suppose elastic channel), the quantum 
mechanical effect of other included channels (through basis set) comes to play 
due to angular momentum coupling through conservation of good quantum 
numbers 'J' and 'M'. 

Two interacting atoms mutually induce a symmetrical pair of dipoles 
[33]. The induced polarization potential is the same for electron-atom and 
positron-atom scattering, in both cases they are attractive [34] . Ps is a highly 
polarizable atom. In coupled channel methodology, the inclusion of the effect 
of dipole polarizability of an atom is possible through excitation of the atom 
to different p-orbitals[10, 15-18, 35-37]. The inclusion of different excited s- 
states provides a short-range strong-coupling force [10] whereas the excited 
p-states supply both a short-range as well as a long-range forces [10,15-17,26- 
27]. The mean static interaction provides a null contribution if Ps. Earlier 
[15] we investigated Ps-H scattering using target-inelastic CCAs to study 
the effect of H(2s) and H(2p) states on elastic channel using target-inelastic 
CCAs. Investigations performed by Burke et al [10,35-36] on e^ — H and 
e^ — H systems established that the inclusion of only the first excited p-state 
of H i.e. H(2p), in the CCA scheme can include 65.77% effects of dipole 
polarizability. Similar properties are expected from Ps as an isotope of H. 
To include the total effect of dipole polarizability one has to choose a CCA 
basis with an infinite number of p-states of both the atoms; it is definitely 
impossible. So one can increase the accuracy of the calculation by expanding 
the basis set only. However, according to the concept of quantum mechanical 
approach, the next target-elastic Ps(3p) channel is most useful. The first 
order effect due to static dipole polarizability vanishes in Ps-atom systems. 
On the otherhand CCA is highly efficient to include the static as well as 
dynamic or non-adiabatic effects [36,38]. Non-adiabatic effect is important 
[17,38-40] if highly polarizabable atom (e.g. Ps). 

The electrons are indistinguishable particles. Two electrons can inter- 
change their positions, the phenomenon is known as exchange. At lower 
energies the effect of exchange is highly important when more than one elec- 
tron in the system. The spins of the target and the projectile electrons can 
make a distribution of 1/4 possibility of forming a singlet (+) state (com- 



bined spin = 0) and 3/4 possibility of forming a triplet (-) state (combined 
spin =1). Accordingly the space part of the system wave function should 
take a symmetric and an antisymmetric form to make the total wave function 
of the system an antisymmetric one. 

We consider a 3-channel projectile-inelastic basis as [H(ls)+Ps(ls,2s,2p)] 
with all the nine possible primary 

Ps{ls) + H{ls) ^ Ps{ls) + H{ls) 
Ps{1s) + hIis) ^Ps(2s) + /J(1s) 
Ps(ls) + nils) -^ Ps{2p) + H{ls) 

and intermediate channels 

Ps{2s) + H{ls) ^ Ps{ls) + H{ls) 

-^ Ps(2s) + H{ls) 

-^ Ps{2p) + H{ls) 
Ps{2p)+H{ls) ^ Pslls) + hIis) 

^ Psl2s) + hIis) 

-^ Ps{2p) + if (Is) 

channels. It is the best 3-channel basis according to concept of quantum 
mechanical approach of eigen-state expansion. All the possible Coulomb in- 
teractions are treated exactly for all the direct and exchange matrix elements. 
Earlier study [26] using such a methodology neglected the most important in- 
termediate channel [Ps(2p)+H(ls)^Ps(2p)+H(ls)] which is responsible for 
non-adiabatic or dynamic effects [38-40]. 

2 Theory 

The total wavefunction of the system \1'^ satisfying the Schrodinger equation: 
i7^±(rp, n, ra) = ^^^(rp, ri, ra) (1) 

is expressed as 

V^ nthuplpLJ^JM ^1 ^^2 Pi 

^^^^ [m m, mJ ( Ml m m) YLMA^i)Yi^rr.M)YkmAh) (2) 
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with 

2 " P 2 " ^ 2 ' ^ ' I Fp I I ri I I r2 I | Fp - ri | | Tp - r2 | ' | ri - r2 
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(3) 

Here P12 stands for the exchange operator and Ri = |(rp + ri) and pi = 
Tp — Fi; i=l,2. ri and r2 are the position vectors of the electrons belonging 
to Ps and H respectively and rp, is that of the positron with respect to the 
center of mass of the system. Unth (r)/r and Vnpip {p)/p are the radial parts of 
the wavefunctions of H and Ps respectively and Fror(k, k', R)/R is the radial 
part of the continuum wave function of the moving Ps atom; Fq indicates all 
the quantities ntltUplpLJiJM of P at the initial channel. 

Projecting the Schrodinger Eqn.(l) just like the Hartree-Fock variational 
approach and integrating over the desired coordinates, we can get a set of 
integro-differential equations. These integro-differential equations can be 
transformed into the integral equations applying the asymptotic boundary 
conditions like Lippmann-Schwinger [14]. These coupled integral equations 
can be formed either in momentum space or in configuration space. Fraser 
et al used the configuration space approach whereas we are using momentum 
space formalism [14]. The set of coupled integral equations obtained for the 
scattering amplitudes are as follows: 

f ± fk' kl - B^ fk' kl - — V [^." ^^'1''^"^"^^'^ ^")fn'n",ni(M", k) 

(4) 
B^ indicate the Born-Oppenheimer [41] scattering amplitudes; plus(+) is for 
the singlet channel for which the space-part of the system wave function is 
symmetric, and minus(-) is for triplet channel for which the space-part of the 
system wave function is antisymmetric. The formulation for Born matrix 
element is available in our previous articles [42-46] and Born-Oppenheimer 
in recent articles [47,48]. Similarly /^ indicate the unknown scattering am- 
plitudes for the singlet and the triplet channels respectively. The summation 
over n"l" is to include various channels. 

These three dimensional coupled integral equations involving fn'i',ni(M-', k) 
and Bn'i',niO^', k) can be reduced to the corresponding one dimensional forms 
through partial wave analysis using the expansion like: 

r^i'n.',nimi^',^) = ikkr'^'J2Il E E E i^'lp, M'^TU^lJ^M^) {Lip, MLmp\J,M,) 

J,M L,Ml L'M'l Ji'^1 J'lM'i 



(J(/;M(m;|JM)(JiZ,,Mimi|JM)Fi,M,,(k')lLM,(k)r-^±(A;'n;n;/;/;,A;npni/p/t) (5) 

Here [A] is the scattering matrix oi NxN type formed by Born and Born- 
Oppenheimer [41] scattering amplitudes, [X] is the column matrix of iV x 1 
type formed by unknown CCA scattering amplitudes and [B] is again a 
column matrix of A x 1 type formed by the Born and Born-Oppenheimer [41] 
scattering amplitudes; the dimension(N) depends on the number of channels 
included in the expansion basis. We have calculated all the Born and Born- 
Oppenheimer amplitudes exactly following an analytic approach and then 
the rest part is carried on following numerical approach using computer and 
FORTRAN programming. The two sets of one dimensional coupled integral 
equations of scattering amplitudes in momentum space for the singlet (+) 
and triplet (-) channels respectively, are solved separately for each partial 
wave(L). 

3 Results and discussion 

We find exact resonances both in singlet and triplet channels in Ps-H scat- 
tering below inelastic threshold in s-wave elastic phase shifts as well as in 
cross sections. The singlet resonances add more importance on the predicted 
binding [1] . The triplet resonance predicts a new binding. Attached Figure 1 
shows our s-wave elastic phase shifts below inelastic threshold and Figure 2, 
the corresponding s-wave partial cross sections. The phase shifts and cross 
sections, both are derived directly from the elastic scattering amplitudes. 
Both the singlet (-I-) and triplet (-) results are presented using the static- 
exchange approximation, the 2-channel projectile-inelastic CCA (CCA2), the 
3-channel projectile-inelastic CCA (CCA3) and using CCA3 but excluding 
the 9th intermediate Ps{2p) + H{ls) -^ Ps{2p) + H{ls). In the first method, 
only the elastic channel is considered in the CCA basis. In the second method, 
the first two primary channels are taken with the intermediates and in the 
third method, all the nine primary and intermediate channels come to play. 
In figure 3, both the singlet and triplet phase shifts of the present CCA3 
are compared with other theoretical investigations of Campbell et al [49], 
Drachman et al [1,50] and Sinha et al [26]. The singlet results using CCA2 
and CCA3 and the triplet results using CCA3 are in complete agreement with 
earlier investigation [26] at lower incident energies near to zero but differing 
near and above the resonance energy region. No phase shift data are available 
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Figure 1: The s-wave elastic phase shifts for Ps-H scattering below threshold. 
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Figure 2: The s-wave elastic cross sections for Ps-H scattering below thresh- 
old. 
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Figure 3: The comparison of s-wave elastic phase shifts for Ps-H scattering 
below threshold. 
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Figure 4: The total cross sections for Ps-H scattering below inelastic thresh- 
old. 
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at this energy region near to threshold for the singlet channel [26,49,51]. The 
triplet result is showing good agreement with other theories but differing 
only near and above the resonance energy region. 

In figure 4, the integrated cross sections using present CCAs are compared 
with the static-exchange and the 22PslH coupled-pseudostate R-matrix re- 
sults of Belfast group. Our integrated/total cross sections are also in good 
agreement with the 22PslH coupled-pseudostate R-matrix results [51] below 
the energy region of triplet resonance. The deviation of the present features 
at the energy region near and above the triplet resonance welcomes a special 
attention. 

The existence of a bound system of Ps and H formed by the electrostatic 
attraction of e"^ and H^i^So) in ground (singlet) state was earlier predicted 
by Drachman [1] and confirmed by others [49-58]. It is first time observed 
using a CCA theory. The present study assures that the cause of this singlet 
resonance is the short-range strong-coupling force arising due to the inclusion 
of Ps(2s) channel with exact exchange and in complete agreement with ear- 
lier prediction. The addition of Ps(2p) channel with exact exchange shifted 
the position of resonance and broadened the width; these findings are also 
consistent with the existing physics. The partial cross section is satisfying 
a formulation like Breit-Wigner which is consistent with a value of non- 
resonant phase shift of ~ 0.5 radian near the resonance energy. The effect of 
this resonance is also reflected on our total/integrated cross sections in figure 
4. The kink in the total cross section curve near the singlet resonance was 
also observed by Sinha el al [26] in their 2-channel CCA calculation which 
is present in our both the CCA2 and CCAS total cross sections in figure 4. 
The present findings of singlet resonance in quite agreement with the reality 
ascertain the accurateness of our present calculation and the computer code. 

The resonance in the triplet (-) channel occur in the CCAS method. It 
signifies clearly the cause of this resonance as the inclusion of Ps(2p) channel 
in the CCA scheme. The inclusion of this channel can provide almost 2/S 
effect [10,S5-S 7] of high polarizability of Ps in a dynamic way [S6]. If we 
exclude the intermediate channel involving target-elastic Ps{2p) -^ Ps{2p) 
transition, the triplet resonance disappears which are displayed in figures 1 & 
2 as 'CCAS ex Ps(2p-2p)'. The s-wave triplet cross section at the resonance 
energy region follows a similar resonance rule for cross section like Breit- 
Wigner and consistent to the value of non-resonant phase shift ~ — 7r/2. The 
effect of this resonance is again reflected in total/integrated cross sections in 
figure 4 as a pronounced dip or a minimum in curve. The importance of non- 
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adiabatic effects in the scattering by a highly polarizable target is discussed 
in Refs. [38-40]. It was estabhshed on comparing the resuhs obtained from 
the non-adiabatic theory by Walters [38] and the CCA theory by Burke and 
Taylor [36] that the CCA theory is quite capable of including the dynamical 
effects by its own nature. 

In addition, the present triplet phase shift data for Ps-H scattering us- 
ing the present 3-channel projectile-inelastic CCA fit nicely [59] with non- 
resonant part as 

Co = -1.4053 + 0.1295E - 0.04613^^ 



and provides the width F = 0.15173eV and resonance position Er = 3.2630eV^. 
A discussion on various existing theories should be useful. A coupled state 
representation of the system wavefunction including the effect of important 
channels is most useful for reliable information. R-matrix theory used by 
Belfast group [49,51] with a coupled pseudostate formalism of the system 
wavefunction is highly useful in this respect. In this theory, the system 
wavefunction is treated accurately upto a certain radius(r < a) and above 
this radius(r > a) the system wavefunction is treated in an approximate way. 
But the effect of the long-range interactions are taken into consideration as 
potentials into the Hamiltonian. The CCA theory used by Calcutta group 
[12-13,15-18,26-27] and the Sao Paulo group [28] is an approach where the 
full configuration space is used through the momentum space formalism [14] 
to get proper long-range effect. The difference should be noted between the 
long-range effect and the long-range interactions; the entire configuration 
space should be used for long-range effect with exact interaction but one 
can approximate the effect of long-range interactions following a trick. The 
convergence problem due to presence of sine and cosine terms in the con- 
figuration space approach to solve the coupled equations at higher angular 
momentum is discussed by Eraser [11,29-30] and Eraser et al[31-32]. The 
convergence is easily obtainable using a momentum space formalism. The 
Sao Paulo group [28] has considered the effect of exchange in an approximate 
way using a model potential. The methodology used by Calcutta group seems 
to be more appropriate; we adapt this methodology. The complex rotation 
and stabilization methods followed by NASA group [1,49, 53-55] and Taiwan 
group [56-57] using Feshbach formalism is a highly efficient tool to detect a 
Feshbach resonance. Other groups of the world e.g. Australia group[60-63]. 
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UCL theory group [64], Italy group [65] are using different variational meth- 
ods. 

An interesting question is the reality of these resonances i.e. if we increase 
the basis set whether the resonances still exist. In atom-atom scattering van 
der Waals interaction is important which comes into picture through excita- 
tion of both the atoms to excited p-states. In principle, the major quantum 
mechanical effect of channel-coupling comes from the closest excitation chan- 
nels i.e. H(ls)+Ps(ls)^H(ls)+Ps(2s) and 

H(ls)-|-Ps(ls)— s>H(ls)+Ps(2p). The target excitation channel and the excita- 
tion of both the atoms channel which is responsible for van der Waals interac- 
tion are energetically far away from the elastic channel H(ls)-|-Ps(ls)^II(ls)+Ps(ls) 
So the effect of these target excitation and both excitation channels can not 
be more important in Ps-H scattering. In addition, the reality of singlet 
resonances are certain. 

4 Conclusion 

To conclude, we introduce a complete treatment of the best 3-channel projectile- 
inelastic CCA for Ps-H scattering considering both the atoms in ground states 
at the initial channel. The scattering phenomena below inelastic threshold 
have been studied. We see resonances in both the singlet and triplet chan- 
nels in the s-wave elastic phase shift as well as in s-wave elastic cross sec- 
tion. Both the resonances are again reflected on the total cross section as 
a peak and a well. Our singlet resonances support the previously predicted 
singlet resonance [1]. The very interesting newly found triplet resonance pre- 
dicts a binding; the cause appears to be a non-adiabatic/dynamical force. 
More investigations using larger basis sets are useful to better understand 
the fact and the mechanism. The minima in s-wave partial cross section for 
triplet channel and in integrated/total cross section may indicate a condensa- 
tion/clustering of the system constitutents. Our theory should be applicable 
to slowly moving atoms and a dilute system in which the inter atomic sepa- 
ration is much greater than the atomic dimension. 
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